Hydrodynamics and the fluctuation theorem 
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The fluctuation theorem is a pivotal result of statistical physics. It quantifies the probability of 
observing fluctuations which are in violation of the second law of thermodynamics. More specifically, 
it quantifies the ratio of the probabilities of observing entropy-producing and entropy-consuming 
fluctuations measured over a finite volume and time span in terms of the rate of entropy production 
in the system, the measurement volume and time. We study the fluctuation theorem in computer 
simulations of planar shear flow. The simulations are performed employing the method of multipar- 
ticle collision dynamics which captures both thermal fluctuations and hydrodynamic interactions. 
The main outcome of our analysis is that the fluctuation theorem is verified at any averaging time 
provided that the measurement volume exhibits a specific dependence on a hydrodynamic time 
scale. 
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The fluctuation theorem states that, in a system 
in steady-state conditions, the probability of observing 
some amount of entropy consumption relative to the 
probability of observing the same amount of entropy 
production is exponentially suppressed. More specifi- 
cally, [IHS] 



P(Slr = -R) 



-Rt 
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where R is the average value taken by the entropy pro- 
duction rate Q over the averaging time = R = 



Among numerous attempts to verify the validity of the 
fluctuation theorem both directly ([T]) and in one of its 
integral forms [4]4TQ], it has been observed that the fluc- 
tuation theorem in its direct form ([T]) does not appear to 
hold in the presence of hydrodynamic interactions [11], 
although, of course, it is still satisfled in an integral form. 

In this contribution, we demonstrate the validity of 
the fluctuation theorem in the presence of hydrodynamic 
interactions, and show that the origin of the previously 
observed apparent inconsistency has a well-deflned phys- 
ical interpretation which lies in correlations and hydro- 
dynamic transport within the fluid. More speciflcally, we 
show that in the presence of hydrodynamics, the entropy 
production rate is a global effect because of correlations. 
When these correlations are eliminated, the extensivity 
of entropy production is lost. Our results not only serve 
as a test of the fundamental laws of thermodynamics, but 
also provide crucial insight for future experiments aimed 
at putting the fluctuation theorem to the test. 

Following Ref. [11 , we consider steady-state planar 
shear flow. We choose shear to be in the x direction with 
a shear rate 7 = dvx/dy. Since in this case the entropy 
production rate can be expressed in terms of the stress 
tensor a^y^ ^ = cr^y^V /ksT ^ where F is a volume and 
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TABLE I: List of setups for which simulations have been per- 
formed together with the respective simulation parameters. 
Setup 7^ is constructed such that hydrodynamic interactions 
(HI) are not present. The viscosities (77) are calculated using 
the analytical formulas of Ref. 13 . For setup 7^ the viscosity 
cannot be evaluated in this way due to the implementation 
of a non-hydrodynamic system. The shear rates are given in 
units of 70 = 0.00625^^^^7^^. 



ksT denotes the thermal energy, the fluctuation theorem 
can be written as [3^ 



V{t,V) = -tV, 
where we have defined the probabihty function 

-A) 
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A key statement of the fluctuation theorem is that 
V{r^ V) is independent of the speciflc properties of the 
system, i.e. the viscosity or the shear rate. Rather, 
7^(r, V) depends on an extensive variable V and a mea- 
suring time r. The extensive variable V is related to the 
volume in which measurement of entropy production is 
performed. 

To study the fluctuation theorem ([2| , we perform com- 
puter simulations employing the method of multiparti- 
cle collision dynamics (MFC). This method is known 
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to correctly capture both thermal fluctuations and hy- 
drodynamic interactions on coarse-grained length- and 
time-scales [121 US] • It has been successfully applied to 
a wide range of steady-state problems, i.e. colloids in 
shear flow [M], polymers in shear flow [15], vesicles in 
shear flow [16] and colloidal rods in shear flow [17]. In 
MPC, the solvent is modeled as a set of point particles 
with continuous coordinates and velocities. The system 
evolves in discrete time steps. At each step, the sol- 
vent particles undergo propagation and collision. During 
propagation, the coordinates fi of each particle are up- 
dated according to its velocity Vi and a time interval dt^ 
U ^ Vi -\- ViSt. During the collision step, all particles 
are sorted into collision cells, and within each cell their 
velocities are rotated with respect to the center of mass 
velocity by an angle a around a random axis. To guaran- 
tee Galilean invariance, the structure of the collision grid 
is shifted randomly before each collision step. The time 
interval dt denotes the time between successive collisions 
within the solvent and is therefore directly related to the 
transport coefficients of the solvent. Thus, we denote 5t 
as the collision time. 

The general framework of calculations is as follows. 
Given a steady-state evaluation of the viscous stress ten- 
sor (Jxy{t)^ where each measurement of the time series is 
performed within a measurement volume V and over a 
measurement time r, we construct the probability distri- 
bution P{cFxy) of all the values in the time series. We 
then average the probability distribution over 100 sim- 
ulation runs with different initial conditions. It is then 
possible to test the ffuctuation theorem 

We fix the simulation parameters to values typical of 
collective, fiuid-like behavior [18 . As normalization, we 
choose the solvent particle mass m = 1, the collision cell 
size a = 1 and the thermal energy ksT = 1. Mass is 
measured in units of m, length in units of a and time 
in units of to = y^ma'^ /ksT. The simulation box sizes 
L = Sa and L = 16a are intentionally chosen small to 
enhance fluctuations. The collision angle a is taken to 
be 130°. The remaining parameters are summarized in 
Table [T] for the several setups for which simulations have 
been performed. 

We employ Lees-Edwards boundary conditions to gen- 
erate planar shear flow . In order to achieve steady- 
state conditions, the system is kept at constant temper- 
ature by a stochastic cell- level thermostat [20 . All simu- 
lations run for 5000 steps before measurements begin to 
be made to guarantee steady-state conditions. 

As the flrst step of the present work, we test an integral 
form of the fluctuation theorem. To this end, we employ 
a simple method to relate the viscosity of the fluid, r^, to 
the rate of temperature increase in the system undergoing 
shear flow in the absence of a thermostat ^21j, 

Cv^=V^' (4) 




FIG. 1: (Color online) Left: thermal energy measured in se- 
tups 2-5 with (t < 5000) and without (t > 5000) a thermostat. 
Right: viscosity measured through the rate of temperature in- 
crease of the fluid (boxes) compared to the exact analytical 
results [13 (solid lines). 



where Cy denotes the volumetric heat capacity. We ob- 
serve that in the absence of a thermostat, the tempera- 
ture of the system increases linearly (Fig. [ij left). The 
results for the viscosity (Fig. [l] right) are in excellent 
agreement with the exact analytical result [T3^ which, in 
our case, gives the shear viscosity as a function of the sol- 
vent density p and the collision time 6t. It is well-known 
that viscosity determined from the long-time limit of the 
viscous stress tensor is in excellent agreement with the 
theory [22 . Here, we have demonstrated that the relation 
between entropy production and a transport coeflicient 
holds, conflrming that overall entropy production is cap- 
tured quantitatively [23j. This result is not a priori evi- 
dent despite the existence of an H- theorem for MPC [24] , 
We remark that the novel method for the determination 
of viscosity using the rate of temperature increase of the 
MPC fluid undergoing shear flow presented here requires, 
in principle, only very short simulation runs, as opposed 
to the determination of the viscosity from the long-time 
limit of the viscous stress tensor [22 . 

Having verifled the integral form of the fluctuation the- 
orem [25], we proceed to test its direct form To this 
end, we measure entropy production by measuring the 
viscous stress tensor a^y as the x-component of momen- 
tum crossing a plane of constant y per unit area and 
per unit time at each simulation step [22 . We choose 
the measurement planes in such a way that on average 
they lie in the middle of each layer of collision cells. In 

the limit of a small mean-free path, <C a, the 

measurement of momentum transfer is therefore limited 
to a single layer of collision cells, since in the collision 
step momentum transfer is conflned to a single collision 
cell. The single-layer measurements of the viscous stress 
tensor, a^^, are then averaged to provide measurements 
over two layers (cr^^), four layers (cr^y) and six layers 
i^xy)' We denote the corresponding probability distribu- 
tions P^^{cF^!y) and the probability functions V^^{r^V). 

The measurement of a^^, where ni is the number of 
measurement layers, corresponds to a measurement of 
entropy production in a volume n^L^a/2. This is due to 
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FIG. 2: (Color online) Top left: measured probability distri- 
butions P{(Txy) for setup 6 for one (green), two (red), four 
(blue) and six (orange) layers over the shortest averaging 
time T — 5t are not zero-centered and not symmetric. They 
become increasingly peaked for increasing measurement vol- 
umes, as expected. Top right: the corresponding probability 
functions P(r, V) normalized by the measurement time t — 5t 
and the volume of a single layer are constant and cor- 

rectly reflect the number of measurement layers. Bottom: 
the resulting scaling of the observation volume V{t) defined 
in Eq. ([5| with the number of measurement planes ni is as 
predicted by the fluctuation theorem for all the simulation 
setups. Different symbols denote different setups. 



the fact that not only is the measured momentum trans- 
fer limited to a single layer of collision cells, but within 
the layer the change of momentum in the two halves sep- 
arated by the measurement plane differs only in sign. We 
test this prediction by computing an observation volume 
V(r) as 



V{t)=V'''{t, V)/t. 



(5) 



In the fluctuation theorem holds, then V^^ (r, V) should 
be constant as a function of a^^. 

We first test these predictions for an averaging time 
corresponding to a single collision time, r = 6t. The 
results are shown on Fig. |2j We find that the probability 
ratio computed from the probability distributions (top 
left) is indeed a constant (top right) and equal to the 
volume of measurement in all cases (bottom). Thus, in 
the short-time limit, the fluctuation theorem holds. 

Next, we compute the running averages of crj^(t) over 
different time frames r G [2St; bOOSt], yielding a function 
^xy,ri'^)^ where r denotes the averaging (measurement) 
time. As before, we express the measurement results in 
terms of the observation volume V{r). 

The results for V{r) obtained from simulations are 
shown on Fig. [s] (top). For the hydrodynamic setups 
(1-6) we find a surprising result. The observation vol- 
ume V{r) is not constant, at least for small values of r. 
This implies that the relation ([2| does not hold in this 



regime. The observation volume scales from the volume 
within which the measurement is performed to the total 
system volume . We investigate this scaling in terms of 
hydrodynamic momentum transport, as this is the only 
transport mechanism in our systems. Hydrodynamic mo- 
mentum transport can be associated with a time scale 
related to the propagation of momentum across a dis- 
tance equal to the lateral simulation box size r/ = L'^ /^^ 
where v is the kinematic viscosity of the solvent. By 
scaling V{t) by the total volume and time by the hy- 
drodynamic time r/, we observe that, for a given number 
of measurement planes, all curves collapse onto the same 
universal curves (Fig.jsj bottom). Therefore, while in the 
short- and the long-time limits the fluctuation theorem 
is satisfied, the volume entering Eq. ([2| is different. In 
the short-time limit it corresponds to the measurement 
volume, whereas at long times it corresponds to the to- 
tal volume of the system. At intermediate times when 
V(t) <C I/^, V{t) exhibits a r^/^ scaling behavior indica- 
tive of simple diffusion in the ^-direction, y'^ ~ t, and the 
fluctuation theorem is also satisfied. 

On the other hand, in the absence of transport phe- 
nomena, one expects the observation volume to be equal 
to the volume of measurement for all averaging times. 
We test this by performing simulations without hydro- 
dynamic interactions. To this end, we choose to adopt 
the following scheme. After each collision step, the ve- 
locities of the particles are randomized. The velocity 
components in the off-shear-directions are sampled from 
a Gaussian distribution with mean zero and variance 
yjk^rjm. The velocity component in the shear direc- 
tion is sampled from a Gaussian distribution with mean 

^Vxl^y {Vi — L/2) and variance ^/kBT/m. Here yi de- 
notes the ^/-component of the coordinate of the solvent 
particle. This procedure guarantees a correct shear gra- 
dient without introducing correlations within the solvent. 

The results for the non-hydrodynamic case, setup 7"*", 
are shown on Fig. [2] and Fig. |3] alongside the results for 
the hydrodynamic setups 1-6. Indeed, we observe that 
while the scaling of the observation volume with the num- 
ber of planes ni is the same (Fig.|2| in all cases, the non- 
hydrodynamic setup does not exhibit any scaling of the 
observation volume with the averaging time (Fig. [3|, in 
sharp contrast to the hydrodynamic case. This is a direct 
consequence of the causality principle: in the absence of 
transport phenomena, it should not be possible to make 
conclusions about parts of the system which lie outside 
of the physical volume of measurement. 

We have shown that the fluctuation theorem ([T]) holds 
in the presence of hydrodynamic interactions at all aver- 
aging times for solvents of different densities and viscosi- 
ties undergoing shear at different shear rates. However, 
transport phenomena have a very significant influence on 
the measurement of fluctuations. In particular, we have 
demonstrated that the interpretation of results pertain- 
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FIG. 3: (Color online) The observation volume V (r) as a function of the averaging time r (top) and the normalized observation 
volume V(t^II? as a function of the normalized averaging time r/r/ (bottom) for different setups, denoted by different symbols. 
The solid lines indicate the total system volumes . The dashed lines indicate the single-, two- and four-layer volumes L^nia/2. 
Shown are measurements of the viscous stress tensor over one plane (left), two planes (middle) and four planes (right). 



ing to measurements of entropy at short times and within 
small volumes in systems with hydrodynamic interactions 
would depend strongly on the properties of the solvent 
and the measurement time. Furthermore, we expect that 
for dynamic probes, i.e. a Brownian particle suspended 
in a solvent, an even finer analysis will be required, as 
here the observation volume will depend not only on the 
properties of the solvent, but also on the properties of 
the probe, such as its diffusion coefficient. We expect 
that, since these results are a direct consequence of the 
causality principle, they are readily generalizable to other 
systems which feature correlations in combination with 
transport phenomena. 
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